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SUMMARY 


A simplified procedure has been developed to calculate the 1*P 
oscillating aerodynamic thrust loads on single -rotation propellers in 
pitch at zero yaw on tractor airplanes. The application of this procedure 
requires only a knowledge of the upflow angles at the horizontal center 
line of the propeller disk. These angles may he obtained by computational 
methods thereby obviating the need for experimental flow surveys. 

The 1-P thrust loads and 1-P vibratory stresses computed by the 
simplified procedure are compared with those calculated using all the 
flow-field characteristics measured for several wing-nacelle -fuselage 
combinations. These comparisons show good agreement for the 1-P thrust 
loads and the 1*P vibratory stresses. 

An analysis is included to show the effects of the various flow-field 
parameters on the 1-P thrust loads. 


INTRODUCTION 


An important consideration in the design of propellers is that which 
deals with vibratory stresses which occur one cycle per revolution (com- 
monly referred to as 1-P) . These stresses are due to 1-P oscillating 
aerodynamic thrust loads imposed on propeller blades as a result of pro- 
peller thrust-axis inclination and/or asymmetries of the flow fields in 
which the propellers operate. Currently, the structural design criterion 
for propeller blades is based upon the allowable amplitudes of the 1«P 
vibratory stresses. Therefore, there is need for a simple, rapid method 
of estimating the amplitude of the 1-P th 2 rust loads. 

It has been shown in previous work (refs. 1 and 2) that steady-state 
propeller theory may be applied to predict adequately the oscillating 
portions of the thrust load on a propeller -blade element around the disk 
if all the characteristics of the flow field in which it is operating 
are known. The amplitude of the 1-P component of this variation (the 1-P 
thrust load) then may be obtained and this load used to determine the 1-P 
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stresses in the blade. This method is quite laborious, requiring com- 
putations for at least l6 points around the disk for each radial position 
to define accurately the variation of the oscillating thrust load. 

Further, the characteristics of the entire flow field at present may be 
obtained only by experimental survey, a cianbersome task which requires 
considerable time. 

This report presents a procedure for calculating the I'P thrust load 
by considering only the upflow angles at the horizontal center line of the 
propeller disk. These angles may be computed theoretically by the methods 
presented in references 3> 5? and thus the need for experimental 

survey is obviated. The procedure is evaluated by comparisons of the 1*P 
thrust loads computed by the procedure with those computed by a more 
rigorous method for several wing-nacelle-fuselage combinations. 


NOTATION 


B number of blades 


c blade chord 


ct 


blade-section thrust coefficient, 


section thrust 
pn^D^ 


ct 


I 


Act, 


blade-section thrust coefficient due to section lift force only, 
section thrust due to l i f^ Appendix A.) 

incremental values of the coefficients, the value at any angular 
position minus the mean value 


D propeller diameter 

K Goldstein correction factor for a finite mmiber of blades 


n propeller rotational speed 

r distance along any radial line from the propeller thrust axis 

R propeller radius 

local velocity at any point (r,fi) at the propeller plane 
(Direction of this velocity is determined by the angles 
Q and See fig. 1.) 
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Vo free-stream velocity 

V’ component of the local velocity in the plane perpendicular to 

a radial line (See fig. 1.) 

— velocity ratio at any point (r,n) at the propeller plane 

''o 

X radial location of any blade section, — 

R 

1*P the fundamental component (one cycle per propeller revolution) 

of the oscillation of the section thrust coefficient or of 
the vibratory blade stresses 

oc^ geometric angle of attack of the airplane thrust axis 

O 0 j_ propeller-induced angle of inflow 

P section blade angle 


7 tan"i [ ^lQ>^^"S^<^tion drag j 

\b lade-section lift/ 

amplitude of the 1-P component of the variation of incremental 
section thrust coefficient due to section lift force, further 
defined as the 1-P thrust load 

6 angle of outflow, measured from a line parallel to the thrust 

axis in a plane through the thrust axis (See fig. 1.) 

0* angle at which the local velocity at any point on the propeller 

disk is inclined to the plane perpendicular to the radial 
line through that point (See fig. 1.) 

p mass density of air in the free stream 

a propeller -blade -element solidity, 

27tr 

^ angle of rotational flow (an apparent, not an actual, rotation) 

measured from a line parallel to the thrust axis in a plane 
perpendicular to a radial line (See fig. 1.) 

(The value of \|r at the horizontal center line of the pro- 
peller disk is identical to the upflow angle at that point 
and corresponds to the notation "A” in ref. 6.) 

angular position about the thrust axis, measured counterclockwise 
from the upper vertical position as seen from the front 
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Subscripts 


\|f effects resulting from consideration of variations in \|r only 

Yi/Yq effects resulting from consideration of variations in 9 and 
Yi/Yo only 

T effects resulting from consideration of variations in \|f, 0, 

and V^/Vo acting simultaneously 


PROCEDURE 

Development 


-Total 


Ct, 


Examinations of total oscillating thrust loads for propellers 
operating in the flow fields of several wing-fuselage-nacelle combina- 
tions have indicated a simple procedure for determining the amplitude of 
the 1*P thrust load. As may be seen in the following sketch of a typi- 
cal total oscillating thrust load and its harmonic components, no 
significant odd-order sinusoidal components other than the fundamental 
(I'P) are present in the variation. i The lack of significant hi^er 

odd-order components in this 
variation (and in all other 
variations examined) makes pos- 
sible the extraction of the 1*P 
component of the thrust load by 
a simple operation which may be 
applied regardless of the ampli- 
tudes or number of even-order 
components present. The sim- 
plicity of this operation stems 
from the fact that all even-order 
harmonics repeat their magnitude 
and sense at half-cycle intervals 
of the fundamental; whereas odd- 
order harmonics, including the 
fundamental, repeat their magni- 
tude but are of opposite sense at 

half -cycle intervals of the fundamental. Thus, if the higher odd -order 
harmonics are nonexistent, one half the arithmetical difference between 
the values of the total oscillating thrust load at any two angular posi- 
tions 180 ° apart (i.e., half -cycle intervals of the fundamental) will yield 
the exact magnitude of the 1-P component at those points. 



90 180 

Angular position, D, 


270 

degrees 


360 


^Harmonic analyses of some of the variations have shown third- and 
fourth-order components of small amplitude but they are of little signif- 
icance in the total variation. 
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The operation may be simplified by utilizing any four equally spaced 
positions around the propeller disk. For the specific points, 0 °, 90 °, 
l80°, and 270°, the expressions shown below may be used to determine the 
amplitude and phase angle of the I’P component. 


1 

2 



J2=900 



f2=2yo°- 


= bc-t^cos cp^ 


( 1 ) 


1 

2 



a=o° 



ii=l 80 ^J 


bct^sin 9 a 


( 2 ) 


where 

bc-t^ the amplitude of the !•? component of thrust load 

9g^ the angle between the position of the maximum 1 *P magnitude and 

the 90° position. 

The small value of 9^^ as shown in the foregoing sketch makes plausible 
the approximation 


Sct^cos 9^ « bc-t^ 

hence, only the value of c^^ for the 90° and 270° positions would be 

required. The procedure then consists of the solution of equation (l) 
rewritten to include the foregoing approximation, as follows: 


ct 


n=900 


- °t 


'I 


fl= 270 °-J 


= Set 


I 


(3) 


While it is necessary to know the values of all of the flow parameters 
at these points to calculate the corresponding values of ct , certain 
characteristics of the flow-field parameters make possible further simpli- 
fications. It is seen in the 
sketch that at the 90° and 270° 


positions, 0 and V^/Vq are 
near 0 and 1.0, respectively, 
while >lf is near its minimum ,^,9 
and maximum values. Thus, the 
values of 0 and Vi/Vq may be 
approximated as 0° and 1.0, 
respectively. Since the values 
of \lf at these positions 
may be computed by methods 
described in references 3 to 6, 
experimental survey of the flow 
field would be unnecessary. 
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It may occur to the reader to circumvent the approximation 

5c^,cos cp^ &c^^ 

hy utilizing expression (2) and the flow-field characteristics at the 0 
and 180° positions. However, it has been found that higher odd-order 
components which are insignificant in the total ^"t^ variation did not 

greatly affect the accuracy of expression (l), but their presence did 
greatly affect the accuracy of expression (2) for the small phase angles 
cp involved. The simultaneous solution of expressions (l) and (2) with 
complete flow-field data did not give greater accuracy than the procedure 
to be examined in this report. 


Evaluation 


As a means of indicating the accuracy and general applicability of 
the procedure, 1>P thrust loads computed by the use of this procedure were 
compared with those determined by a l6-point solution. Measured flow-field 
data were used for the l6-point solution. ^ For the procedure, measured 
values of t were used and the values of 0 and V^/Vq were assumed to 
be 0° and 1.0, respectively. The blade stresses (due to aerodynamic load- 
ing alone) resulting from the 1-P thrust loads computed by the two methods 
were also compared. The geometries of the model configurations considered 
in these comparisons are shown in figure 2. For purposes of this evalua- 
tion, it was assumed that a four-bladed right-hand propeller was installed 
on each model. The characteristics of this propeller and the method 
employed in the computation of the stresses are given in reference 1. 

The 1*P thrust loads for the 0.7-radial station of the propeller 
blade for all models as computed by the two methods are compared in fig- 
ure 3. Comparisons of the radial variations of the l-P thrust loads 
computed by both methods for two models are shown in figure 4. The I’P 
blade stresses for these two models are shown in figure 5. From these 
figures, it is seen that the procedure is generally satisfactory. 


DISCUSSION 


Although the procedure of this report may be considered sufficiently 
accurate for the configurations investigated, it was concluded that it 


^The values of the flow parameters used in these computations were 
determined by experimental survey and were reported in references 2, 3, 
and 6. Since data were obtained at the port propeller planes, the flow 
fields for the starboard propellers were assumed to be the reflected 
images of those for the port propeller. 
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would be impossible from a cursory examination of a configuration to 
predict whether this accuracy could always be realized. To gain an 
insight into the cause or causes of the 1*P thrust loads, it was decided 
to investigate the separate effects of the various flow-field parameters. 
Therefore, the variations of the total thrust loads on a blade element 
were calculated for the models in the following manner: (l) using meas- 

ured values of f, and assuming 0=0° and V^/Vo = 1.0, (2) using meas- 
ured values of 0 and V 2 /V 0 and assuming \|f = 0 , (3) using measured 
values of f, 6, and From these variations, the 1*P components 

were extracted. As may be seen from figure 1, the first case shows the 
effect of the angular variation of the inflow velocity V» , while the 
second case shows the effect of the magnitude variation of the inflow 
velocity V» . While it is not apparent from the equations for the total 
thrust variations (see Appendix A) that this method would yield the 
exact effects of the variables, it was found that the summation of the 
thrust loads for cases (l) and (2) was nearly identical with the thrust 
variation calculated for case (3). The 1*P components of thrust varia- 
tions and the associated flow-field parameters are shown in figure 6. 

From these figures, it is seen that \|r produced a 1-P component 
which peaked at approximately 90° and 270°. Since the procedure of this 
report utilizes only f at the 90° and 270° positions, it essentially 
determines the 1*P thrust load due to \|/. From these figures, it is also 
seen that 0 and V^/Vq produce a significant 1*P component which peaked at 
approximately 0° and l80°. The models which have a large l‘P component due 
to 0 and V 2 /V 0 are characterized by a relatively large blunt nacelle. 
Further, when the nacelle is mounted on the wing, the magnitudes of these 
1-P components are increased due to the velocity gradient induced by the 
^ing. It should be noted, however, that all of these 1*P components due 
to 0 and V^/Vo are so phased as to add little to the total 1-P com- 
ponent. The total 1-P component appears to be more dependent on the phase 
angle between the 1-P components due to \|r and 0 and Vi/Vq than on the 
magnitude of the 1-P component due to 0 and V^/Vq. 

Although the foregoing analysis has shown the causes of the 1-P 
thrust load, the inability to predict the phase angles involved makes it 
impractical to use this analysis to attempt refinement of the basic pro- 
cedure of this report. The analysis may prove useful for estimation of 
probable accuracies of the basic procedure for other configurations. 


CONCLUDING REMARKS 


A procedure has been presented for calculating the 1-P thrust load 
on propellers of tractor airplanes in pitch at zero yaw, whose applica- 
tion requires only a knowledge of the upflow angles at the horizontal 
center line of the propeller disk. The evaluation of the procedure 
showed that the 1-P thrust loads and the 1-P blade stresses resulting 


8 


NACA TN 3395 


therefrom could be predicted for the airplane configurations investigated 
and it is believed satisfactory results could be obtained for similar 
configurations . 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif,, Mar. 21, 1955 
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APPENDIX A 


The 1*P stress problem is primarily flatwise bending of the blade 
which results from oscillations in the blade normal force. Therefore^ it 
is the sum of the components of the thrust and torque forces perpendicular 
to the chord line of the blade element which must be determined. The 
normal forces on a blade element may be expressed as 

N -|3 = t cos p + fqsin p (A1) 


where 

t section thrust 

fq section torque force 

The section thrust on a blade element was given in coefficient form in 
reference 2 as 


c^ = 




57.3 / . 

(^cot cp + 


V^sin 
jrnDx j 


(A2) 


where 


cp = cpo + o-i 


cp„ = tan-/ 

V^imDx - V’ sin i)f ) 

This expression may be simplified to consider only the section thrust 
due to lift. 


"tz 


Kjt^x^ 


0-1 cot cp 

^ “if 


sin \lf 
TrnDx 


(A3) 


This is the expression which was used to calculate the values of c^ 
considering all the flow-field parameters as measured values. The values 
of 0° and 1.0 may be substituted for Q and Vz/Vq in this expression as 
described in the text. The resulting expression is 


10 


NACA TN 3395 


c-t-j = Kit^x^ 


cot cp 


0-1 

57-3 , 

cot Cp + ^7^ 


1 - 


VoSin \lr 
itnDx 


where 


cp = cpo + ai 


cp = tan'll 


VqCOb y 


omDx - Vo sin \|c 


If the blade normal force is expressed in coefficient form and both sec 
tion thrust and section torque force are expressed in terms of section 
thrust due to lift, cp, and y , equation (Al) becomes 


■=Nb = 8 + ct,(l - |^)tan(9 + r)sln S 

When equation (4) is rearranged 


cos 3 


tj fcos[(y - cp) - 2(P - cp)] + cos (7 + cp) 


cos (7 - cp) + cos (7 + cp) 


(A4) 


(A5) 


For simplicity, the term in the brackets was assumed to be unity which 
is tantamount to saying that the drag force on a blade element does not 
contribute materially to the blade-element normal force. Therefore, 
only the thrust due to lift on the blade element was considered to obtain 
the normal force which is expressed as 

r. 

Nb ” cos P 
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A-14137.2 

Figure 1.- Schematic representation of the flow field parameters. 
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Nacelle coordinates 

Station 

Radius 
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21.5 

/ 

23.0 
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24.0 
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25.0 
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26.0 

9 

27.0 
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26.5 

24 
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36 

29.6 

42 

29.7 

46 

29.8 

54 

29.6 

62 

29.8 

72 

29.8 

96 

29.5 

120 

29.0 

144 

26.5 

168 

22.4 

192 

16.4 

217 

n .2 


AH dimensions are in inches. 


(a) Model a 


Figure 2.- Geometric characteristics of the models. 
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(b)Mode! b 
Figure 2 — Continued. 


16 


NACA TN 3395 


AH dimensions are in inches 



Wing characteristics 

Aspect ratio 9.9 

Section normal to c/4 line 64A0I0 

Area 265 ft ^ 

Sweep of c/4 line 
Wing incidence 
Taper ratio 
Scale 


y33.0^ 



307.2 


(c) Model c 

Figure 2.— Continued. 
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A// dimensions are in inches 



Aspect ratio 

73 

Section normal to c/4 line 

64 A 010 

Area 

233 ft 

Sweep of c/4 tine 

40^ 

Wing incidence 

3.4 ** 

Taper ratio 

0.49 

Scale 

0.4! 


Propeller disk 


64.8 



247.2 


76.8 



[40Y / 



/ 

/ 

/ 

/ 


-1- 

o/ 

I 

^ tOiC.U 




i 

384.0 


(d) Model d 


Figure 2. — Concluded. 


IP load due to ^ only at Sl=90^and 270°, 
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O 1 ' ^ ' 

0 .04 .08 .12 .16 


hP load due to the entire flow field, (SCf^)j 

Figure 3. — Comparison of IP loads as computed by the procedure of 
this report with those computed using the entire flow field for 
several wing-fuselage-nacelle configurations; x-0.7, Vo=!65 mph, 
n- 1250 r pm, /S^.^ff2l.7° 



Fraction of tip radius ^ x 
(a) Model b, port propeller. 

Figure 4. — Comparison of the radial variation of the IP load computed by the 
procedure of this report with that computed using the entire flow field; 
Vo=l65mph, n^!250rpm, /Sqj‘21.7^, = 10^. 
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(b) Model b, starboard propeller. 
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fcj Mode! c, port inboard propeller. 
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• 3 .4 .5 .6 .7 .8 .9 t.O 


Fraction of tip radius, x 
(a) Model b 

Figure 5 — Comparison of IP vibratory stresses resulting 
from I P loads computed by the procedure of this report 
with those resulting from IP loads computed using the 
entire flow fields Vo=l65mph, n=!250rpm, ^^^^=21.7;^ aQ=IO? 


IP vibratory stress, psi 
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(b) Model c 
Figure 5— Concluded. 


Angle of outflow, S, 
Angle of rotational flow, 
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Angular position, a, deg 
(a) Models a and b 

Figure 6. — The flow- field parameters and their etfects on the TP variation of incremental section thrust 
coefficient for several wing- fuselage— nacelle combinations; Vo=!65mph, n= !250rpm, 
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(b) Mode! c, inboard 
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Figure 6— Continued. 
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(c) Model c, outboard 
Figure 6. — Continued. 
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Figure 6. — Concluded. 


